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ABSTRACT 

We present a comparison between the SCUBA Half Degree Extragalactic Survey 
(SHADES) at 450 /zm and 850 /zm in the Lockman Hole East with a deep Spitzer 
Space Telescope survey at 3.6 — 24 /im conducted in Guaranteed Time. Using stacking 
" analyses we demonstrate a striking correspondence between the galaxies contributing 

I^S , the submm extragalactic background light, with those likely to dominate the back- 

' grounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photomet- 

ric redshifts, we show that at least a third of the Spitzer-identified submm galaxies 
at 1 < z < 1.5 appear to reside in overdensities when the density field is smoothed 
at 0.5 — 2Mpc comoving diameters, supporting the high-rcdshift reversal of the local 
star formation - galaxy density relation. We derive the dust-shrouded cosmic star 
formation history of galaxies as a function of assembled stellar masses. For model 
stellar masses < lO^^M©, this peaks at lower redshifts than the ostensible z ^ 2.2 
maximum for submm point sources, adding to the growing consensus for "downsiz- 
ing" in star formation. Our surveys are also consistent with "downsizing" in mass 
assembly. Both the mean star formation rates {dM^,/dt) and specific star formation 
rates {{1 / M^,)dM^, / dt) are in striking disagreement with some semi-analytic predic- 
tions from the Millenium simulation. The discrepancy could either be resolved with a 
top-heavy initial mass function, or a significant component of the submm flux heated 
by the interstellar radiation field. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: star-burst - infrared: galaxies - submillimetre 



1 INTRODUCTION 



The SCUBA Half Degree Extragalactic Survey (SHADES, 
Mortier at al. 2005, Coppin et al. 2006) is a long-term submm 
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Figure 1. The SHADES 850 /xm map with point sources sub- 
tracted, used for stacking analyses in this paper (greyscale). Also 
plotted is the coverage of the Spitzer 24 fim data (solid red line) , 
the 3.6 fim data (solid blue line), 4.5 fj,m (solid green line), 5.8 
(dashed red line) , and 8 (dashed magenta line) . Note the avail- 
ability of Spitzer data over most of the SHADES area. 
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Figure 2. Correlation between model stellar mass and absolute K 
magnitude predicted using the Millenium simulation by de Lucia 
(2006), for several redshift shells. Also plotted are the data from 
the MUNICS survey (Drory et al. 2004). Note the small dispersion 
and the clear redshift dependence. 



survey conducted at the James Clerk Maxwell Telescope 
from 2003-2005. A key goal of SHADES has been to de- 
termine whether submm galaxies are the likely progenitors 
of giant ellipticals. The clustering of submm galaxies is a 
strong discriminant of competing models (van Kampen et 
al. 2005), and measurements of the angular correlation func- 
tion of submm galaxies in broad Az ~ 0.5 redshift shells is 
one of the principal experimental aims of SHADES. In this 
paper we will take a different approach to the problem, by 
estimating the matter overdensities in which submm galax- 
ies reside via the assembled stellar masses in the submm 
galaxy environments. The SHADES survey was conducted 
in two fields, each with abundant multi-wavelength support- 
ing survey data. The ~ 0.1 deg^ surveyed by SHADES in 
the Lockman Hole East field, in particular, has some of the 
best Spitzer Space Telescope data of any contiguous field 
over hundreds of square arcminutes. The comparison be- 
tween SHADES and this Spitzer data, which was taken in 
Spitzer guaranteed time, forms the basis of our constraints 
on the submm galaxy environments, and allows us impor- 
tant new insights on the submm extragalactic background 
light. 

The galaxies that dominate the extragalactic back- 
ground light at any given redshift are necessarily the same 
as those which dominate the comoving volume- averaged lu- 
minosity density at that redshift (e.g. Peacock 1999). The 
favourable K-corrections in the submm make the submm 
extragalactic background light sensitive to the cosmic star 
formation history throughout most of the history of the Uni- 
verse. Resolved submm point sources from blank field sur- 
veys (i.e. those at the few-mjy level) contribute a few tens of 



z 0.025 1.05 2.05 3.05 4.1 5.25 

a 2.007 1.457 1.156 0.8193 0.3064 0.5831 
b -0.3944 -0.4010 -0.4058 -0.4178 -0.4373 -0.4190 

Table 1. Best-fit values for the conversion between K-band AB 
absolute magnitude and stellar mass using equation [T] from the 
de Lucia (2005) Millenium simulation. The top row gives the red- 
shift, and the other rows give the fitted parameters in a narrow 
redshift bin centred on that redshift. 

percent to the 850 fim extragalactic background, but cannot 
account for all of it (Hughes et al. 1998, Barger et al. 1998, 
1999, Blain et al. 1999, Eales et al. 2000, Scott et al. 2002, 
Smail et al. 2002, Cowie et al. 2002, Scott et al. 2006). At 
350 — 450 /im there are very few reliably detected resolved 
point sources (e.g. Scott et al. 2002, Khan et al. 2005, 2007), 
and those that have been detected are far from account- 
ing for the majority of the 450 nm extragalactic background 
light. However, 850 /.im-selected galaxies can be readily de- 
tected at 350 — 450 fim and can account for a minority of the 
350 — 450 /im background (e.g. Chapman et al. 2005, Khan 
et al. 2005, 2007, Kovacs et al. 2006, Coppin et al. 2007). 

This situation, particularly at 450 pm, will change with 
the advent of the SCUBA-2 camera on the James Clerk 
Maxwell Telescope. In the meantime, attention has focussed 
on stacking analyses. Instead of aiming to detect individual 
resolved galaxies, this approach seeks to detect the average 
signal from a population. This has met with some success. 
Peacock et al. 2000 found a ~ 3cr signal at 850 /im from the 
Lyman-break population in the Hubble Deep Field North. 
The submm:UV fiux ratio suggested an obscuration very 
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Figure 3. Model stellar mass estimates of Spitzer galaxies dis- 
cussed in the text, versus photometric redshifts from Dye et al. 
(2008, SHADES paper VII). Note that the model stellar masses 
above z = 0.5 depend only weakly on photometric redshift. 
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Figure 4. Density distribution percentiles for submm galaxies 
(large symbols), compared with that of 3.6/im-selected galaxies 
(dots). The matter density distribution estimated from Spitzer 

3.6 /irn imaging in broad redshift bins (indicated by dashed lines) 
was smoothed with a 0.5 comoving Mpc diameter top hat kernel 
as discussed in the text. Regions with zero estimated density were 
excluded from this test. Note that submm galaxies sample the 
highest 10 — 20 percentiles of this matter density distribution at 
lower redshifts. 
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Figure 5. Density distribution percentiles for submm galaxies, 
compared with that of 3.6 ^m-selected galaxies (dots). The mat- 
ter density distribution estimated from Spitzer 3.6 /im imaging 
in broad redshift bins (indicated by dashed lines) was smoothed 
with an 2 comoving Mpc diameter top hat kernel as discussed in 
the text. Regions with zero estimated density were excluded from 
this test. Note that around half of the submm galaxies in the lower 
redshift bin sample the highest 10 — 20 percentiles of this mat- 
ter density distribution, when smoothed to this scale. The larger 
smoothing width reduces the area with no density field data, so 
the sample of submm galaxies is larger. 



different to that of the submm point source population, and 
there were hints of a flat redshift distribution at z > 1 in 
these faint submm-emitting galaxies. Further stacking anal- 
yses of extremely red galaxies (e.g. Webb et al. 2004, Takagi 
ct al. 2007) found them to contribute a significant minority 
of the obscured star formation history. A submm stacking 
analysis of near-infrared and mid-infrared selected galaxies 
from the 5' x 5' Spitzer Early Release Observations (Serjeant 
et al. 2004) found that Spitzer 5.8 /im and 8/^m populations 
could account for around a (juartcr of the 850 /nm extra- 
galactic background light, and the majority of the 450 /jm 
background, albeit in a small sample. The sample size was 
not large enough in this study to distinguish the stacked 
signal from low-rcdshift red dusty galaxies, and that from 
high-rcdsliift galaxies. 

To constrain the redshift ranges responsible for the 
submm extragalactic background light, larger samples were 
needed. Wang et al. (2006) and Dye et al. (2006) both made 
stacking analyses of Spitzer-selccted galaxies, though in the 
former case it was combined with H-band selection. Both 
groups found that Spitzer galaxies contribute significantly 
to the submm extragalactic background light. However, the 
redshift ranges responsible in these surveys differed, with 
Wang et al. finding the z < 1 population dominating (their 
figure 12), while Dye et al. found z > 1 the more important 
(their figure 6) though with slightly larger errors. One pos- 
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sible explanation for this difference is cosmic variance; the 
wide-area SHADES would be ideal to resolve this contro- 
versy. Another possibility is the effect of redshifted PAH fea- 
tures in the Dye et al. analysis, which would not be present 
in the < 3.6/im Wang et al. analysis. 

In this paper we extend these results to a deeper 
Spitzer catalogue, and a wider-area submm survey. In a 
confusion-limited submm survey, the stacking signal-to- 
noise is roughly proportional to the square root of the num- 
ber of submm beams, and since SHADES is the widest-area 
contiguous submm survey to a depth approaching the effec- 
tive point source extraction limit (Scott et al. 2002), this is 
the best opportunity to date to examine the submm stacking 
signal of Spitzer galaxies. 

This paper is the ninth in the SHADES series of papers. 
Paper I (Mortier et al. 2005) presented the survey design, 
motivation and data analysis. Paper II (Coppin et al. 2006) 
presented further data analysis, the source counts, the cat- 
alogues and the maps. Paper III (Ivison et al. 2007) gave 
the radio and Spitzer 24 /xm identifications of the submm 
galaxies in SHADES. Paper IV (Aretxaga et al. 2007) made 
photometric redshift estimates of the SHADES catalogue 
galaxies using the far-infrared to radio spectral energy distri- 
butions. Paper V (Takagi et al. 2007) examined the submm 
properties of near-infrared galaxies in the SHADES suvey 
data in the Subaru-XMM Deep Field. Paper VI (Coppin 
et al. 2007) presented the results of 350 fj,m observations 
of a subset of SHADES sources. Paper VII (Dye et al. 
2008) made fits to the spectral energy distributions of the 
SHADES galaxies in the Lockman Hole, and paper VIII 
(Clements et al. 2007) performed a similar analysis for the 
SHADES galaxies in the Subaru-XMM Deep Field. Paper 
X (van Kampen et al. 2007) measures the clustering of the 
submm galaxies in the SHADES survey. A further series of 
papers will concern the 1.1mm data taken to supplement 
the SHADES survey with the AzTEC instrument on the 
James Clerk Maxwell Telescope. 

This paper is structured as follows. The Spitzer and 
submm data are summarized briefly in section [2] Section [3] 
describes our methodology and results. We discuss the con- 
text of our results in section |4j and we draw conclusions 
in section [5] Throughout the paper we assume a "concor- 
dance" cosmology, with density parameters Qm — 0.3 and 
Qa ~ 0.7, and a Hubble constant of Ho = 72 km s~^ Mpc"^. 



2 DATA ACQUISITION 

The submm data were taken from 2003-2005 at the James 
Clerk MaxweU Telescope (JCMT) with the SCUBA cam- 
era, in submm opacities of 0.265 < rgso^m < 0.283 and 
1.41 < T450Mm < 1.52, i.e. JCMT weather bands 2-3. Chop- 
ping/nodding was performed at position angles of and 90 
degrees, with chop throws of 30", 44" and 68", though at 
the centre of the map only the 30" /90 degree combination 
was used (Scott et al. 2002). The submm data acquisition, 
calibration, reduction and analysis are described in full in 
Mortier et al. (2006). A noise- weighted point source filtering 
was made on the maps, and the separate chop/nod images 
were combined optimally (Serjeant et al. 2003a, Mortier et 
al. 2006). Note that the SCUBA 450 /im absolute flux cal- 
ibration is typically uncertain to ~ 30%. The maps from 



this analysis are used in this paper. Three further parallel 
data reduction efforts are described in Coppin et al. (2006) , 
and 850 fim maps from some of these reductions have been 
used to test the robustness of the results presented in this 
paper. The 450 /xm maps used in this paper are from the 
SHADES data reduction described as the "primary" reduc- 
tion in Mortier et al. (2005), and analysis "B" in the later 
Coppin et al. (2006). These maps have the lowest noise of 
the available 450 fim maps, and have short timescale opacity 
variations modelled using the water vapour meter. No reli- 
able 450 ^m point sources are detected in SHADES, though 
this is not to say that one cannot find point source candi- 
dates in the maps. A cross-comparison of the four reduc- 
tions found very few overlaps between the candidate source 
lists; this demonstrated that the 450 /im opacity during the 
SHADES runs (see appendix A of Coppin et al. 2006), and 
the stability of the opacity, were not suited to reliable point 
source extraction. We will show that this does not preclude 
statistical constraints on the 450 /im-emitting populations. 
The 850 /xm maps of the four data reduction methods are in 
excellent agreement (Coppin et al. 2006), and the SHADES 
point source list is derived from the consensus of the four 
analyses. 

The Spitzer data were taken in Guaranteed Time, using 
the IRAC and MIPS instruments (Fazio et al. 2004, Rieke 
et al. 2004). As shown in figure [fl most of the SHADES 
survey area in the Lockman Hole East is covered at 3.6 /im 
(4.47/.iJy, 3cr), 4.5 /.im (4.54/iJy, 3cr), 5.8 /xm (20.9MJy 3cr), 
S fim (12.5/jJy, 3.2o-) and 24/im (38/xJy, 4a). We only use 
Spitzer IRAC sources detected in at least two Spitzer bands. 
The area was also mapped at 70 /^m and 160 /im, and the 
comparison between this data and SHADES will be the sub- 
ject of a future paper (Egami et al. in preparation). Most 
of the SHADES Lockman field was also covered by a deep 
15 fim survey with the CAM instrument on the Infrared 
Space Observatory, further details of which can be found 
in Elbaz et al. 1999 and Rodighiero et al. 2004. We select 
galaxies with 15 fxm flux densities above 100/^Jy for this pa- 
per. The Subaru-XMM Deep Field was also observed by 
SHADES, and has Spitzer data from the SWIRE survey 
(Lonsdale et al. 2004), but since this data is significantly 
shallower than the Lockman data we do not consider it here. 

BRIz imaging was obtained from the SUPRIMECAM 
instrument on the Subaru telescope, to 5a point source 
depths of 26.8, 25.8, 25.7 and 25.0 in B, R, I and z respec- 
tively (3" diameter AB magnitudes). K-band photometry 
was obtained from the UKIRT Deep Infrared Sky Survey 
(UKIDSS, Lawrence et al. 2007) to a point source sensitiv- 
ity of 22.9 {5a AB magnitude). Further details are in Dye 
et al. (2008, paper VII). 



3 METHODOLODY AND RESULTS 

3.1 Photometric redshift estimates 

We use the Spitzer galaxy photometric redshift catalogue 
of Dye et al. (2008). This catalogue is derived using the 
HYPER-Z code (BolzoneUa et al. 2000) applied to the nine- 
band optical-IRAC photometric catalogue. Further details 
of the spectral energy distribution templates are found in 
Dye et al. (2008) . In figure |3] we show the model stellar 
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Figure 6. Histograms of 850 /im flux in the SHADES Loclcman map as a whole (unhatched bars) compared to the 850 /xm fluxes at 
the positions of the Spitzer/ISO galaxies selected at the indidcated wavelengths (hatched bars). The left-hand vertical scales show the 
numbers for the hatched histograms, and the right-hand scales refer to the unhatched histograms. Note that in all plots, the hatched 
bars lie typically to the right of the unhatched bars. The results of a comparison using the Kolmogorov-Smirnoff test are also shown in 
the diagrams. Submm point sources have been removed. 



masses (derived below) as a function of the photometric 
redshifts. Note that redshift aliasing can scatter galaxies to 
erroneously high redshifts (and hence to high masses), and 
there are plausible examples of this in figure O We have 
opted not to impose any arbitrary cuts in the photomet- 
ric redshift catalogue to remove these outliers, and instead 
leave this to the discretion of the reader. Redshift aliasing 
and consequent erroneously high stellar masses in a small 
subset would not alter the statistical conclusions of this pa- 
per. 

Dye et al. (2008) also uses optical and Spitzer pho- 
tometery to derive photometric redshift estimates for the 
SHADES galaxies themselves, and comparisons with other 
redshift estimators can be found in Dye et al. (2008), Aretx- 
aga et al. (2007) and Clements et al. (2007). We wiU use the 
Dye et al. (2008) determinations in this paper. The main dis- 
advantage of the photometric redshifts in the Lockman Hole 
is the lack of spectroscopic training sets, though a compari- 
son of the Dye et al. photometric redshifts of submm-selected 
galaxies with their spectroscopic redshifts, and between in- 
dependent photometric redshift determinations (their fig- 
ures 3 and 4) shows the photometric redshifts are accurate 



to |A2:|/(1 + 2:) ~ 0.09 consistent with other studies (e.g. 
Chapman et al. 2005, Pope et al. 2006). This is more than 
sufficient for our purposes. 



3.2 Mass estimates of Spitzer galcixies 

The 3.6 /^m and 4.5 /im IRAC bands are dominated by the 
redshifted hght from old stellar populations, and are there- 
fore useful estimators of the assembled stellar masses. The 
K-corrections are also simple at these wavelengths, since 
it samples the Rayleigh- Jeans tail of the stellar photo- 
spheric emission. We use the model spectral energy distri- 
butions from Dye et al. (2008) to obtain rest-frame K-band 
monochromatic luminosities; the results are insensitive to 
the assumed spectral energy distribution. 

There is no accepted conversion between rest-frame K- 
band luminosity and stellar mass, as a function of red- 
shift. Our aim is to compare our results with the de Lu- 
cia (2005) Millenium simulation, and one can only do this 
self-consistently by adopting a conversion consistent with 
that simulation. Therefore, we adopt an empirical conver- 
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Figure 7. Histograms of 450 /im flux in the SHADES Loclcman map as a whole (unhatched bars) compared to the 450 fj,m fluxes at 
the positions of the Spitzer/ISO galaxies selected at the indidcated wavelengths (hatched bars). The left-hand vertical scales show the 
numbers for the hatched histograms, and the right-hand scales refer to the unhatched histograms. Note that in all plots, the hatched 
bars lie typically to the right of the unhatched bars. The results of a comparison using the Kolmogorov-Smirnoff test are also shown in 
the diagrams. 



sion based on these simulations, based on fits to the simu- 
lated data in figure [S] This conversion is consistent with the 
observed luminosity-dependent evolution in K-band stellar 
mass-to-light ratios in the MUNICS survey (Drory et al. 
2004). Our conversion is 

logio M* = a{z) + b{z) X /Cabs.AB (1) 

where the values of a and b are interpolated from the best 
fit values tablulated in table [T] 

In calculating the matter overdensities, we add a dark 
matter contribution following the total mass estimates from 
gravitational lenses by Ferreras et al. 2005, who found that 
the total mass scales as the stellar mass to the 1.2 power, 
with total mass equalling stellar mass at 3.18 x IO^Mq. Be- 
low this mass we make no correction for dark matter con- 
tribution. Our results are not sensitive to the dark matter 
assumptions. 

3.3 The density field around submm galaxies 

As shown in Blake et al. (2006), the SHADES survey does 
not have enough sources or field galaxy redshifts to accu- 



rately determine the galaxy-SHADES cross-correlation func- 
tion. We therefore used an alternative estimator which is es- 
sentially a stack of the 3.6 ^m Spitzer data at the positions of 
the SHADES sources. Using Spitzer galaxy mass estimates 
discussed above, we created projected mass density maps in 
broad redshift bins, and smoothed each map with a top-hat 
circular kernel with 0.5 — 2 comoving Mpc diameters at the 
central redshift of the bin. We then compared the projected 
matter density at the positions of the SHADES galaxies with 
the histogram of matter density for the map as a whole, us- 
ing methods similar to established techniques for stacking 
analyses (see e.g. section \3A^ . We omitted SHADES galax- 
ies lying in regions with zero density (no 3.6 fim identifica- 
tion nor sufficiently close neighbours) and also restricted the 
comparison to the non-zero density regions of the map as a 
whole. 

This comparison is shown in figures U and [S] where the 
percentile of the submm galaxies' environment is plotted 
against redshift. In the 0.5 Mpc smoothing case, the submm 
galaxies in the lower redshift bin lie in the top 10 — 20 per- 
centile of the matter density distribution. However, the small 
smoothing kernel leaves large regions of the map with no 
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Figure 8. Signal-to-noise images of stacked 850 /im postage stamps centred on Spitzer/ISO galaxies, as discussed in the text. The 
greyscale is from —3cr to +4(7. The scales marked are arcseconds. Note the clear detections at the centres of all the images. Negative 
sidelobes are also visible. 



density field data, so the sample of submm galaxies is small. 
In the 2 Mpc smoothing case, the sample size is more than 
doubled; only five galaxies are excluded in the lowest red- 
shift bin. The submm galaxies appear to lie in a wide range 
of environments, but around a third of the submm galaxies 
at 1 < z < 1.5 lie in the top 15 percentile of the density 
distribution. 



We performed a similar calculation for 3.6 /xm-selected 
galaxies, which are plotted as small dots in figures |4] and 
[H These appear to have a more uniform distribution in 
the lower redshift bin than the submm galaxies. We com- 
pared the submm galaxy and 3.6 /xm-selected populations in 
the 1 < z < 1.5 bin using a Kolmogorov-Smirnoff test. In 
the 2 Mpc smoothing case the difference is only marginally 
significant (32% probability that the distributions are the 
same), but in the 0.5 Mpc case the distributions are clearly 
different (0.2% probability that the distributions are the 
same). Note that in making this comparison, we have sub- 
tly changed the question we are asking. We wish to know 
if submm galaxies lie in richer environments than average, 
but this average can be taken as a volume average or as a 
per-galaxy average. The raw percentiles address the former, 
and the submm-Spitzer comparision addresses the latter. 



3.4 The submm emission of 3 — 24 Spitzer 
galaxies 

In Serjeant et al. (2004), the regions of the map near submm 
point sources were simply excluded from the analysis. This 
runs the risk of removing the submm signal from compan- 
ions to the submm point source. Ifere, point sources detected 
in the combined point-source-filtered 850 /im map with sig- 
nificance levels of > 3.5cr were subtracted from the orig- 
inal chop/nod images, which were then filtered with the 
chopped/nodded point source kernels, and the maps were 
optimally combined to create a residual 850 fim image. This 
is the same as the procedure adopted by Dye et al. (2006), 
and this is the threshold used for the submm point source 
catalogue in Coppin et al. (2006). (We will show below in 
figures [6] and [7] that our results are not sensitive to this 
threshold.) There is evidence that the submm point source 
population has different Spitzer:submm flux ratios than the 
Spitzer-selected population (Serjeant et al. 2004), so it is im- 
portant to remove the point source population before stack- 
ing. No reliable 450 /im point sources are detected, so the 
point-source-filtered maps at 450 fim are used without mod- 
ification. Our methodology differs from that of Wang et al. 
(2006), in which submm point sources were left in the map. 

The submm point spread function sums to exactly zero, 
because of the negative sidelobes from the chopping and 
nodding. Therefore, there is no risk of overestimating the 
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Figure 9. Signal-to-noise images of stacked 450 /im postage stamps centred on Spitzer/ISO galaxies, as discussed in the text. The 
greyscale is from — Scr to +4(7. The scales marked are arcseconds. Note the clear detections at the centres of all the images. Negative 
sidelobes are also visible. 



submm flux of a given Spitzer galaxy by also counting its 
neighbours (Peacock et al. 2000, Serjeant et al. 2004). This 
is because the expectation value of the submm flux from 
(unclustered) neighbouring galaxies equals that of the map, 
which is exactly zero by virtue of the zero-sum point spread 
function. The effect of the clustering of the Spitzer popula- 
tion is estimated below in section |3l6l 



We can therefore calculate the mean submm flux of 
Spitzer galaxies by averaging the measurements at the po- 
sitions of Spitzer galaxies in the submm maps, even if there 
is > 1 Spitzer galaxy per submm beam. The difficulty in 
extracting 450 /im point sources (section [2| raises the pos- 
sibility of non-Gaussian features in the maps contributing 
to the signal. Such non-Gaussian features could be caused 
by, for example, imperfect sky subtraction or imperfect cor- 
rection for atmospheric opacity; the sky is many orders of 
magnitude brighter than the extragalactic signal (e.g. Ser- 
jeant et al. 2003a). However, by the Central Limit Theorem, 
the probability distribution of the mean of a sample is ap- 
proximately Gaussian with a variance cr^/n where is the 
variance of the underlying distribution being sampled; the 
distribution of the mean is increasingly Gaussian for larger 
samples and for more Gaussian-like underlying distributions. 
Here, the underlying distribution is well-approximated as 
Gaussians (e.g. Mortier et al. 2005), and many thousands 
of samplings are taken from the underlying distribution. We 



are therefore confident that the mean flux levels of stacked 
populations are Gaussian-distributed. 

The mean flux level has the advantage of being phys- 
ically interpretable, but it is not necessarily the most effi- 
cient detection statistic of a stacked signal. The Kolmogorov- 
Smirnoff test has been widely used to test whether the dis- 
tribution of submm fluxes at the positions of interest are 
a random sampling from the map as a whole (e.g. Ser- 
jeant et al. 2004, Dye et al. 2006). This test is asymptot- 
ically distribution-free, and is therefore insensitive to non- 
Gaussian features in the underlying maps. Furthermore, 
a comparison with a control sample (such as randomized 
submm source positions) is intrinsic to the test, since it com- 
pares the flux distribution of the map as a whole with the 
map fluxes at the positions of interest. It is also possible 
to translate the Kolmogorov-Smirnoff significance level into 
an equivalent number of cr of a Gaussian distribution, by 
inverting Pks = 0.5erf(a/V2) + 0.5. 

Some authors have advocated the use of the error- 
weighted mean fluxes, rather than the unweighted mean 
stacked fluxes (e.g. Dye et al. 2006). In Serjeant et al. (2004) 
it was argued that it was not obvious that the Central Limit 
Theorem would apply to these weighted quantities. How- 
ever, Dye et al. (2006) tested their weighted means and 
found them to be only subtly biased. The large area of the 
SHADES maps may provide an advantage: provided the field 
of view is sufficiently large, and provided the noise level is 
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sufficiently uniform, the zero-sum point spread function of 
each individual source will still produce a zero net contri- 
bution on average to the S/N (signal-to-noise) and S/N^ 
images, as well as to the flux image S. In this paper we will 
make error-weighted coadded submm postage stamps of the 
sources to be stacked. We will show that our error-weighted 
mean fluxes show no evidence for a systematic shift relative 
to the unweighted mean fluxes. 

We applied the methodology of Serjeant et al. (2004) 
to test whether the submm fluxes at the positions of 
Spitzer/ISO galaxies were representative of the submm map 
as a whole, or whether there is on average positive flux 
at the Spitzer/ISO galaxy positions. We excluded regions 
of the 850 fim map with noise levels above 5mJy, and re- 
gions of the 450 fim map with noise levels above 20mJy, 
since these disproportionately affect the non-noise-weighted 
stacks. At both wavelengths, we then calculated the median 
noise levels in the unmasked regions, and then masked all 
areas with more than twice these noise levels. The 450 /im 
data quality is much more dependent on the weather condi- 
tions, so more of the 450 pm map is excluded by our noise 
cuts (e.g. table [Sjl . The results are shown in flgures |6] and 
[T] We detect the Spitzer galaxies at 3> 99% confidence at 
both 450 fim and 850 fim, and for all Spitzer wavelengths 
from 3 — 24 fim. We also detect the ISO 15 /xm-selected 
population at 99.6% confidence at 850 /^m and 96% con- 
fidence at 450 ^m. The mean fiuxes are given in table (2] 
as are the Gaussian-equvalent a values of the Kolmogorov- 
Smirnoff significance levels (though recall the uncertainties 
in the 450 fim fiux calibration noted above) . Our stacked 
mean fiuxes are in good agreement with the previous de- 
terminations of Serjeant et al. (2004), though in addition 
we have made clear submm stacking detections at 3.6 
and 4.5 fim. The mean fiux ratios for the 5.8 — 8 ^m-selected 
galaxies presented in tableware somewhat lower than those 
presented in Serjeant et al. (2004) , which may be due to the 
lack of brighter Spitzer sources in the very small field of view 
of the Early Release Observations. There are some difficul- 
ties in interpreting these fiux ratios, as they are summed 
from galaxies spanning a range of redshifts, and the con- 
tribution made by individual galaxies will depend on their 
location in the luminosity-redshift plane. Furthermore, the 
mean mid-infrared flux is sensitive to the presence of the 
few brighter sources in the sample, which may lead to un- 
derestimates in the quoted errors in the flux ratios. We will 
return to this topic in section |4l Another anomalous flux 
ratio is the 450/im:15 /xm ratio; we believe this is due partly 
to small number statistics in this sample, and partly to the 
fact that only brighter mid-infrared galaxies are detected 
at this wavelength which may bias the sample to submm- 
weak AGN dust tori. These galaxes are also only marginally 
detected (table (2]). 

The Kolmogorov- Smirnoff test reports an appar- 
ently unrealistically small significance when comparing the 
850 /xm fiuxes of 24 /im-selected galaxies with the submm 
map as a whole. Examination of the numbers in the bins of 
the histograms shows why this is the case; in effect, the prob- 
ability that the distributions are identical is immeasurably 
small. To our knowledge, this is the best submm stacking 
detection ever made. 

We constructed a noise weighted sum of the postage 
stamps around Spitzer/ISO galaxies selected at each wave- 
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Figure 10. Evolution in comoving volume- averaged star forma- 
tion rates as a function of total mass for 3.6 /^m-selected galaxies, 
in ±0.5 dex mass bins centred on the masses shown. The brighter 
submm point source population discovered by Smail et al. 1997, 
Hughes et al. 1998 and Barger et al. 1998 are not included, be- 
cause we removed submm point sources prior to our stacking anal- 
ysis. Symbols: the circle is 10^ Mq, upward triangles 10^*^ Mq, 
downward triangles IO'^^Mq, diamonds IO^^Mq, and the square 
is the 10^^ Mq constriant. Note that the lower-mass galaxies have 
star formation histories peaking at much lower redshifts than the 
putative 2; ~ 2.2 peak in the bright submm point source popula- 
tion. 



length. Figures [S] and [5] show the signal-to-noise images of 
the galaxies selected at these wavelengths. We clearly have 
strong detections of our sample at all wavelengths. Fur- 
thermore, the off-centre positions in these stacked postage 
stamps provide a control, and confirm the stacking is not 
prone to false positives. The weighted means are quoted in 
table [2] and are in good agreement with the unweighted 
means. In particular, there is no evidence in table [2] for 
weighted means systematically offset from the unweighted 
mean fluxes; any systematic offset must be far smaller than 
the random noise in the measurements. 

Since our signal-to-noise is so high (unlike in the much 
smaller sample of Serjeant et al. 2004), we can investigate 
the sub-populations which dominate the stacking signal. In 
figure [To] we plot the comoving volume-averaged star forma- 
tion rate estimated from the 850 /im stacked fluxes, assum- 
ing an M82 spectral energy distribution shape, and using a 
conversion derived for this SED from the Kennicutt (1998) 



conversion: 



SFR 



(2) 



Me/year 5.8 x IO^Lq 3.6 x IO^Lq 

This assumes a Salpeter initial mass function from 0.1 
to 100 AIq . Our estimator for the total submm flux contri- 
bution from a population of galaxies in a logarithmic mass 
interval AlogjQ A4 is 

1 



Ftot(mJy deg ^dex ^) 



AA logio M 



Er=iM(i)/c(i) (3) 
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Table 2. Mean submm fluxes for Spitzer-selected and ISO-selected populations. Also quoted are the numbers of galaxies in each stacking 
analysis. Kolmogorov-Smirnoff signiflcance values are quoted in flgures |6] and [T] and translated in this table into an equivalent number 
of (T of a Gaussian distribution as discussed in the text, except for the 850 fim detection of 24 fim sources which has a significance level 
too high for numerical inversion. The quoted errors on the mean fluxes are a/\/~N, where A'^ is the number of measurements and the 
variance of the measurements, except in the case of flux ratios where the a/\/N errors in each quantity have been propagated. Known 
submm sources have been subtracted from the maps prior to this analysis, as discussed in the text. The effect of the clustering of Spitzer 
galaxies is discussed in section [3.61 
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Figure 11. Left: mean star formation rates, dM/dt, of our galaxies. Symbols as in flgure lTol Right: corresponding predictions from the 
de Lucia et al. (2006) simulation. Note that the observed star formation rates clearly exceed the predictions. 



where A is the survey area, Ai (i) is the submm map flux at 
the position of galaxy i (of which there are A^) and c(i) is 
the completeness of the Spitzer catalogue for galaxies similar 
to i. The calculation is not sensitive to the completeness 
correction. To first order both the submm and Spitzer fiuxes 
are constant over the redshift intervals in question (due to 
the negative K-corrections at both wavelength ranges) . Our 
estimate of the error on the total flux is 



:V^tl{M{i)/c{i)r 



(4) 



For conversion from 850 /im fluxes to luminosities, we use 
the central redshift of the bin, noting that the submm lu- 
minosity is roughly independent of redshift over the redshift 
ranges considered. We do not plot bins in which galaxies 
with masses equal to the mass at the centre of the bin are not 
above the Spitzer flux limit throughout the bin. We weight 
the flux contributions of each galaxy according to its accessi- 



ble comoving volume. Figure [TT] shows the mean (quiescent) 
star formation rate per galaxy. Figure [T^] shows the mean 
(quiescent) star formation rate per unit galaxy mass, as a 
function of mass. This quantity has the dimensions of one 
over time (e.g. Gyr~^). This characteristic star formation 
timescale can be regarded as the timescale over which the 
bulk of the galaxy's baryonic matter would be converted 
into stars (though for consistency with elsewhere in this pa- 
per we use the total mass estimates for the galaxies, not 
just baryonic). The quantities plotted in figure \TI\ have no 
dependence on the completeness of the Spitzer catalogue, 
though the errors on the quantities depend on the sample 
size. We have also plotted the mass-doubling timescale as 
a function of redshift on this figure, assuming SFR/M, is 
constant; galaxies above this line may be regarded as star- 
bursting. A further useful metric, also plotted in this figure, 
is the specific star formation required to build up the entire 
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Figure 12. Left: star formation rate per unit galaxy mass, (1/ M)dM / dt. This has dimensions of [T]~^ and the corresponding star 
formation timescales are given in the right-hand ordinate. Symbols as in figure 1101 Also plotted is the specific star formation rate 
required to double the stellar mass by z = assuming a constant specific star formation rate, and the rate required to assemble the 
galaxy since the Big Bang assuming a constant star formation rate. Points lying above either line may be regarded as starbursting. Note 
that most of our detections can be regarded as starbursting. Right: corresponding predictions from the de Lucia et al. (2006) simulation. 
Note the obvious discrepancies with the observations. 



observed stellar mass since the Big Bang, assuming SFR is 
constant. For comparison, the specific star formation rates 
of the SHADES point sources are shown in figure [TS] We 
will discusss these figures in section \A72\ 

3.5 The cosmic near-infrared and submm 
backgrounds 

In figures [14] and [15] we plot the contributions to the cos- 
mic submm background Hght made by Spitzer/ISO galax- 
ies, as a function of their near/mid-infrared fiux. We correct 
the Spitzer/ISO catalogues for incompleteness by compar- 
ison with published source counts from Rodighiero et al. 
(2004), Papovich et al. (2004), and Fazio et al. (2004). The 
figures show the submm background contribution per decade 
of near/mid-infrared flux, and compare these contributions 
to those made by the same galaxies to the near/mid-infrared 
backgrounds calculated from the published source counts. 

There is a remarkably strong correspondence between 
the 24 /^m-selected galaxy contribution to the 850 /^m extra- 
galactic background, and to the 24 ^m-selected background. 
There is also a correspondence between the 3.6 — 4.5 /im- 
selected galaxy contributions to the 450 /im extragalactic 
backgrounds, and to the 3.6 — 4.5 /im backgrounds. The 
24 fim background contributions also correlate well with the 
450 /im contributions. We will discuss the reasons for these 
correspondences in section [4] 

In figures [TSl and [T71 we integrate the data in figures [HI 
and 1 151 and plot the cumulative contributions to the cosmic 
submm background hght, as a function of Spitzer/ISO flux. 
It is clear that about a quarter of the extragalactic 850 /im 
background hght is resolved by Spitzer, and the majority of 
the 450 /^m extragalactic background is resolved. 
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Figure 13. Star formation rate per unit galaxy mass, 
(l/M)dM/dt, for the SHADES point sources in the Lock- 
man Hole. As with figure 1121 the corresponding star formation 
timescales are given in the right-hand ordinate. Objects with 
spectroscopic redshifts are shown in filled symbols, and objects 
with optical-Spitzer-based photometric redshifts from Dye et al. 
(2008) are shown as open symbols. Note that the star formation 
timescales for SHADES sources arc much shorter than in the pop- 
ulation as a whole. 
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Figure 14. Contributions of tiic Spitzcr/ISO galaxies to tlic 850 fini cxtragalactic baclcground light (data points, left-hand ordinates) as 
a function of Spitzer/ISO flux. Also plotted are the contributions of the same galaxies to the cosmic near-infrared and mid-infrared back- 
grounds (curves, right-hand ordinates), estimated from the source counts as discussed in the text. Note in particular the correspondence 
between the 24 /im and 850 /xrci extragalactic backgrounds. 



3.6 The impact of clustering on stacking analyses 

Because of the limited resolution provided by SCUBA, we 
should be clear that the signal detected in the stacking anal- 
ysis represents a contribution from any object within around 
10" of the target Spitzer/ISO galaxies. Because these galax- 



ies will be surrounded by a population of correlated neigh- 
bours, it is therefore possible that the stacked flux gives an 
overestimate of the emission from the target galaxies. This 
is simple enough to estimate: the additional flux is just the 
integral of the background intensity, /, times the angular 
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Figure 15. Contributions of the Spitzcr/ISO galaxies to the 450 /im extragalactic background Hght (data points, left-hand ordinates) as 
a function of Spitzcr/ISO flux. Also plotted are the contributions of the same galaxies to the cosmic near-infrared and mid-infrared back- 
grounds (curves, right-hand ordinates), estimated from the source counts as discussed in the text. Note in particular the correspondence 
between the 3.6 — 4.5 //m and 450 //m extragalactic backgrounds. 



cross-correlation between the target galaxies and the back- For a Gaussian beam, and assuming w = {0/9o) this 
ground, w{6), times the beam B: gives S = 0.20(6'o/arcsec)°-*mJy at 850 Hm and S = 

0.28 (0o/arcsec)°'* mJy at 450 nm. The appropriate value of 
00 is of course open to debate, but Oliver et al. (2004) mea- 
P sure ^0 — 1.24" for the 3.6 /im-scloctod population. The 

S = I w{9) B{9)27Td d9 (5) stacked fluxes at 450 /^m are at least five times the max- 



14 Serjeant et al. 



imum value that could arise from neighbours (taking the 
extreme case in which all the 850 /xm signal arises in this 
way). At 850 /xm the contribution may be more significant, 
but the contribution estimated from equation [5] is necessar- 
ily an over-estimate because not all the 850 /xm background 
is attributable to these Spitzer galaxies. We can use figure 
1161 in which only ~ 25% of the 850 ^m background is at- 
tributable to the Spitzer galaxies, to estimate iteratively the 
correlated 850 /xm fiux from the other Spitzer galaxies. This 
reduces the clustered contribution by a factor of four, yield- 
ing 5* = 0.06 mjy with = 1.24", which is at most a 
20 — 30% correction to the 850 /xm fiuxes quoted in table 
(2] While non-zero, this is not sufficient to affect our conclu- 
sions. 



4 DISCUSSION 

4.1 The link between the near-IR and submm 
backgrounds 

Any stacking analysis is only capable of determining the 
first moment of the distribution; the mean values in figures 
[To] and [12] may belie a large variation in the population. We 
have also subtracted point sources, so strongly starburst- 
ing galaxies are omitted from these figures; we are therefore 
probing only the mean quiescent levels of star formation in 
these galaxies. 

Our stacking detections have much higher signal-to- 
noise than any previously obtained, partly because of the 
depths of the Spitzer and SCUBA surveys, and partly also 
because the stacking signal-to-noise scales with the square 
root of the number of submm beams and SHADES has the 
largest contiguous submm survey fields to date. Our 450 /xm 
stacking results are the best indicators of the populations 
that will be found to dominate the 450 jim background by 
SCUBA-2. The prospects are good for follow-ups of the ul- 
tradeep SCUBA-2 Cosmology Survey, because the Spitzer 
galaxies that appear to dominate the 450 yLxn background 
are less challenging targets for 8 — lOm-class spectroscopy 
than SCUBA point sources (e.g. Serjeant et al. in prepara- 
tion). Similarly, the prospects appear good for spectroscopic 
follow-ups of ALMA point sources below the SCUBA-2 con- 
fusion limit. The clustering of bright submm point sources 
as a function of redshift is a key goal of the SHADES survey 
(e.g. van Kampen et al. 2005, Mortier et al. 2006), providing 
strong constraints on semi-analytic models of galaxy evolu- 
tion; similarly, the redshift-dependent clustering of galaxies 
a factor of ~ 10 fainter in submm flux, for which optical 
follow-up is easier, is likely also to provide a strong con- 
straint on such models. 

Paradoxically, it is the brighter submm point sources 
that are the most challenging to follow up. The populations 
sampled by the SPIRE instrument on Herschel will be chal- 
lenging to follow up in the optical (see also e.g. Khan et 
al. 2005, 2007). Spectroscopic redshifts for such populations 
may be better determined in the medium term by molecular 
line spectroscopy (e.g. Wagg et al. 2007), and in the longer 
term by SPICA (e.g. Nakagawa 2004) 

We find that the submm:24^m ratios for most 24 /im- 
selected galaxies are very different to those of most submm- 
selected galaxies, in agreement with Serjeant et al. (2004). If 



the 24 /^m population had submm:24/im fiux ratios consis- 
tent with those of the submm point source population, the 
24 /im sources considered in this paper would overproduce 
the 850 /xm background by a factor of ~ x3. We argue that 
the bulk of the Spitzer population has a quiescent star for- 
mation level much lower than that of submm point sources, 
while the latter are heavily obscured objects (e.g. Serjeant et 
al. 2003b, Clements et al. 2004, Small et al. 2004) that are 
challenging for optical follow-ups, and far-infrared "loud" 
with high specific star formation rates (figure [13)) and short 
star formation timescales. Such episodic star formation is 
supported by models of AGN feedback in massive galaxies 
in the early Universe, and by the small inferred mass accre- 
tion rates onto central supermassive black holes in submm 
point sources (e.g. Alexander et al. 2005). 

The physical sizes of populations dominating the far- 
IR background are rather smaller than those inferred for 
submm point sources. At a fiux density of ~ 15/iJy, the 
8/im population which we have found to contribute signif- 
icantly to the 450 /xm extragalactic background has optical 
identifications in our imaging with typical diameters H 1" 
(the optical identification diameters of submm point sources 
can be up to ~ 2-4", e.g. Small et al. 2004, Pope et al. 2005, 
though other authors claim sub-arcsecond sizes, e.g. Chap- 
man et al. 2004, Biggs & Ivison 2008). The regions of star 
formation in these galaxies will be resolvable with ALMA, 
which will probe the cool large-grained dust phase, and Dar- 
win direct imaging which will probe the transiently-heated 
small grains and PAH phases. This also suggests that 0.1" 
is the coarsest resolution that would be useful for a future 
~ 50 — 200 /im far-infrared interferometer (FIRI) to resolve 
the internal structure of individual galaxies that comprise 
the cosmic far-infrared background. 

4.2 The mass dependence of star formation 

When plotting the total contribution to the submm back- 
ground from redshift shells in our Spitzer samples (not 
shown), we find that the z < 1.5 population is dominant, 
similar to the results of Wang et al. (2006). However, this 
neglects the fact that different luminosity and mass ranges 
are sampled at different redshifts. This may be one under- 
lying cause of the difference between the Wang et al. (2006) 
and Dye et al. (2006) stacking results, since the K-correction 
effects are different in their respective samples; cosmic vari- 
ance is another possibility. Ours is the first direct attempt 
to segregate the mass contributions to the submm-derived 
cosmic star formation history. 

The mass segregation in figure [10] shows evidence for 
star formation in galaxies with model stellar masses ~ 
IO^'^Mq assembling the bulk of their stellar masses at much 
lower redshifts than the estimated z ~ 2.2 peak in the 
submm point source population. This implies an increasing 
dominance at higher redshift of higher mass systems in the 
volume-averaged star-formation rate. These observations are 
in accordance with qualitative expectations from "downsiz- 
ing" in star formation (Matteucci 1994, Bressan et al. 1996, 
Cowie et al. 1996; see also papers VII and VIII, Dye et al. 
2008 and Clements et al. 2007 respectively). 

Figure [TT] shows the mass- and redshift-dependence of 
the mean star formation rates per galaxy, rather than per 
unit volume. The decrease in the mean star formation rates 
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per galaxy at ~ 10^ Mq and 2 < 1 can only be recon- 
ciled with their increasing volume- averaged contribution if 
their number densities are increasing, in agreement with ex- 
pectations from mass downsizing (e.g. Pozzetti et al. 2007). 
Figure [11] also shows the corresponding predictions from the 
de Lucia (2005) Millenium simulation. Our star formation 
rates are a factor of a few higher than those predicted in 
this simulation. 

The star formation timescales in figure [T2] at redshifts 
2 > 1 scale approximately inversely with the mass of the 
system, quite unlike e.g. a Schmidt law for local late-type 
galaxies. In figure [T^] we again compare this to the pre- 
dictions from the de Lucia (2005) Millenium simulations. 
Again, there is a striking discrepancy with the simulation 
predictions. 

How bad would the photometric redshifts have to be in 
order to explain the disagreement between model and data? 
The results in figure [T2] are surprisingly insensitive to pho- 
tometric redshift, and it is therefore not likely that errors 
in the photometric redshifts are the cause of the disagree- 
ment. This insensitivity to redshift errors is due to the model 
stellar mass and star formation rates both being relatively 
insensitive to redshift (e.g. figure |3} because both are to 
some degree subject to negative K-corrections. The specific 
star formation rate estimates therefore depend mainly on the 
submm;near-IR fiux ratio, and not on redshift. It is there- 
fore hard to see how a redistribution of galaxies among the 
redshift bins could bring the data into agreement with the 
models; removing galaxies from one redshift bin to make it 
agree better with the models would make the disagreement 
worse in the other bins into which the galaxies are moved. 

One hint that the star formation rates may be overesti- 
mated comes from asking what will happen to the 10^^ Mq 
galaxies if they continue forming stars at these high con- 
tinuous rates. The galaxies appear to be forming stars at 
over ten times the rate required to double their mass by the 
present day, almost regardless of redshift at 2 < 2. This is 
clearly not sustainable; it would require a rapid truncation 
of the star formation at slightly higher masses, and a con- 
tinual feeding of lower mass galaxies into the ~ 10^" A/0 bin, 
for which there is no evidence in the stellar mass functions 
of galaxies (e.g. Bell et al. 2004, 2007, Caputi et al. 2006). 

There are at least two possibilities that might reduce 
our estimates of the number of stars forming in these galax- 
ies. One approach explored successfully by the Baugh et al. 
(2005) application of the Millenium simulation is to assume a 
top-heavy initial mass function in star-forming galaxies. The 
reduction in the number of stars formed in this model is sim- 
ilar to the discrepancy between the de Lucia (2005) predic- 
tions and our measurements. However, the top-heavy initial 
mass function is normally only applied to extreme starburst- 
ing systems, and not to the quiescent star formation level 
in the galaxy population as a whole. A second possibility 
is that the observed-frame submm fluxes are dominated not 
by star formation, but by cool cirrus heated by the galaxies' 
interstellar radiation flelds. The galaxies would then have 
a cooler SED than the M82 template assumed above. We 
demonstrate the strength of this effect in figure [181 Such a 
model was proposed for bright submm point sources by Ef- 
stathiou & Rowan-Robinson 2003 (see also Clements et al. 
2007, paper VIII), and although mm- wave and radio inter- 
ferometry have not on the whole yielded the large angular 
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Figure 16. Cumumative contributions to the 850 fim extragalac- 
tic background light, as a function of near/mid-infrared flux, for 
various Spitzer/ISO surveys. The horizontal dashed line is the 
total background derived by Lagache et al. (1999), the uncertain- 
ties on which are less than a factor of 2. Note that the Spitzer 
contributions converge to roughly a quarter of the total 850 
background. 
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Figure 17. Cumumative contributions to the 450 fira extragalac- 
tic background light, as a function of near/mid-infrared flux, for 
various Spitzer/ISO surveys. Symbols as figure 1161 Note that 
the Spitzer contributions are capable of accounting for the en- 
tire 450 background, within the errors. 

sizes predicted by these models (e.g. Tacconi et al. 2006, Ivi- 
son et al. 2007), it remains possible that these models are 
broadly correct descriptions of the fainter submm popula- 
tion. 

4.3 The environment dependence of star 
formation 

We have attempted to measure the matter environments of 
submm galaxies. Semi-analytic models predict that submm 
galaxies should be strongly clustered and lie in some of 
the largest overdensities at their redshifts (e.g. van Kam- 
pen et al. 2005) , for which we have found tentative evidence 
(e.g. figures [4] and [5} and which agrees at least qualitatively 
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Figure 18. Bolometric luminosity as a function of redshift for 
an observed 850 flux of 1 mjy, for three spectral energy dis- 
tributions: M82 (full line), a cirrus-dominated spectrum from Ef- 
stathiou et al. 2000 (short dashed, green), and an Arp 220 model 
also from Efstathiou et al. 2000 (long dashed, green). 



with previous measurements (e.g. Blain et al. 2004, Blake 
et al. 2006). Recently, Elbaz et al. (2007) and Cooper et 
al. (2007) have shown evidence in the GOODS and DEEP2 
surveys that z ~ 1 star-forming galaxies (as evidenced by 
24/im emission) are preferentially found in richer galaxy en- 
vironments. This trend is in the opposite sense to that seen 
in the local Universe. Furthermore, the observed environ- 
ment dependence of star formation is stronger than that 
predicted by semi-analytic models. Our observations extend 
this trend to higher redshifts and star formation rates. The 
properties of the star formation density field, as opposed to 
the galaxy density field, may be a key arena for the future 
confrontation of data and semi-analytic predictions. Future 
Herschel and SCUBA-2 surveys will yield large catalogues 
of bright submm-selected galaxies, and their near-infrared 
and submm environments will be easily measurable with the 
warm AKARI/Spitzer missions and SCUBA-2 follow-ups re- 
spectively. 



5 CONCLUSIONS 

There is a strong correspondence between the galaxies that 
dominate the submm extragalactic background light and 
those that are detected in deep Spitzer surveys. The submm- 
derived specific star formation rates in the Spitzer popula- 
tions are much higher than those predicted by some semi- 
analytic simulations; this may be due to a component of 
submm emission heated by the interstellar radiation fields 
leading to overestimates of the star formation rates, or to a 
top-heavy initial mass function in the Spitzer galaxies, or to 
some unknown deficiency in the models. We find evidence for 
downsizing in both star formation and mass assembly. We 
also find evidence that around a third of submm-selected 
galaxies at redshifts 1 < z < 1.5 lie in the upper ~ 20%-ile 
of the galaxy density distribution, in contrast to the redshift 
zero tendency of star-forming galaxies to avoid the richest 
environments. 



ACKNOWLEDGEMENTS 

We would like to thank the referee, Wei-Hao Wang, for a 
careful reading of the manuscript and for helpful comments. 
SS would like to thank the Science and Technology Facil- 
ities Council for support under grants PP/D002400/1 and 
PP/D003083/1. KC and SD would like to thank the Sci- 
ence and Technology Facilities Council for support. IRS ac- 
knowledges support from the Royal Society. AP acknowl- 
edges support provided by NASA through the Spitzer Space 
Telescope Fellowship Program, through a contract issued 
by the Jet Propulsion Laboratory, California Institute of 
Technology under a contract with NASA. Data was ob- 
tained at the James Clerk Maxwell Telescope under program 
M/02B/U52. The James Clerk Maxwell Telescope is oper- 
ated by The Joint Astronomy Centre on behalf of the Sci- 
ence and Technology Facilities Council of the United King- 
dom, the Netherlands Organisation for Scientific Research, 
and the National Research Council of Canada. This work is 
based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology under a contract 
with NASA. 



REFERENCES 

Alexander, D., et al., 2005, ApJ, 632, 736 

Aretxaga, I., et al., 2007, MNRAS submitted ( |astra-ph/0207503| l 
[Paper IV] 

Barger, A.J., et al., 1998, Nature, 394, 248 

Barger, A. J., Cowie, L.L., Sanders, D.B., ApJ, 518, L5 

Baugh, CM., et al., 2005, MNRAS, 356, 1191 

Bell, E.F., et al., 2004, ApJ, 608, 752 

Bell, E.F., Zeng, X.Z., Papovich, C, Borch, A., Wolf, C, Meisen- 

heimer, K., 2007, ApJ, 663, 834 
Biggs, A.D., & Ivison, R.J., 2008 MNRAS in press 

tarXiv:0712 '.3047 ) 
Blain, A.W., Kneib, J.-P., Ivison, R.J., Small, I., 1999, ApJ, 512, 

L87 

Blain, A.W., Chapman, S.C., Small, I., Ivison R., 2004, ApJ, 611, 
725 

Blake, C, Pope, A., Scott, D., Mobasher, B., 2006, MNRAS, 368, 
732 

Bolzonella, M., Miralles, J.-M., Pello, R., 2000, A&A, 363, 476 
Bressan, A., Chiosi, C, Tantalo, R., 1996, A&A, 311, 425 
Caputi, K.I., et al., 2006, ApJ, 637, 727 

Chapman, S.C., Small, I., Windhorst, R., Muxlow, T., Ivison, 

R.J., 2004, ApJ, 611, 732 
Chapman, S.C., Blain, A.W., Small, I., Ivison, R.J., 2005, ApJ, 

622, 772 

Clements, D., et al., 2004, MNRAS, 351, 447 
Clements, D., et al., 2007, MNRAS submitted [Paper VIII] 
Cooper, M.C., et al., 2007, MNRAS submitted (arXiv/0706.4089) 
Coppin, K, et al., 2006, MNRAS, 372, 1621 [Paper II] 
Coppin, K., et al., 2007, MNRAS submitted [Paper VI] 
Cowie, L.L., Songaila, A., Hu, E.M., Cohen, J.G., 1996, AJ, 112, 
839 

Cowie, L.L., Barger, A.J., Kneib, J.-P., 2002, AJ, 123, 2197 

de Lucia, C, Springel, V., White, S.D.M., Croton, D., Kauff- 

mann, G., 2006, MNRAS, 366, 499 
Drory, N., et al., 2004, ApJ, 608, 742 
Dye, S., et al., 2006, ApJ, 644, 769 

Dye, S., et al., 2008, MNRAS submitted ( ar Xiv:0802.0497l l [Paper 
VII] 

Eales, S.A., et al., 2000, AJ, 120, 2244 



SHADES paper IX 17 

Efstathiou, A., & Rowan- Robinson, M., 2003, MNRAS, 343, 322 
Efstathiou, A., & Rowan- Robinson, M., Siebenmorgen, R,., 2000, 

MNRAS, 313, 734 
Elbaz, D., et al., 1999, A&A, 351, L37 
Elbaz, D., et al., 2007, A&A, 468, 33 
Fazio, G.G., et al., 2004, ApJS, 154, 10 

Ferreras, I., Saha, P., Williams, L.L.R., 2005, MNRAS, 623, L5 
Hughes, D.H., et al., 1998, Nature 394, 241 

Ivison, R., et al., 2007, MNRAS submitted (:,astro-ph/0702544| | 

[Paper III] 
Kennicutt, R.C., ApJ, 498, 541 
Khan, S., et al., 2005, ApJ, 631, L9 
Khan, S., et al., 2007, ApJ, 665, 973 

Kovacs, A., Chapman, S.C., Dowell, CD., Blain, A.W., Ivison, 

R.J., Small, I., Phillips, T.G., 2006, ApJ, 650, 592 
Lagache, G., Abergel, A., Boulanger, F., Desert, F.X., Puget, J.- 

L., 1999, A&A, 344, 322 
Lawrence, A., et al., 2007, MNRAS, 379, 1599 
Lonsdale, C., et al., 2004, ApJS, 154, 54 
Matteucci, F., 1994, A&A, 288, 57 
Mortier, A., et al., 2005, MNRAS, 363, 563 [Paper I] 
Nakagawa, T., 2004, AdSpR, 34, 645 
Oliver, S., et al., 2004, ApJS, 154, 30 
Papovich, C., et al., 2004, ApJS, 154, 70 

Peacock, J. A., 1999, Cosmological Physics, CUP, Cambridge 

Peacock, J. A., et al. 2000, MNRAS, 318, 535 

Pope, A., et al., 2005, MNRAS, 358, 149 

Pope, A., et al., 2006, MNRAS, 370, 1185 

Pozzetti, L., et al., 2007, A&A, submitted {arXiv/0704.1600) 

Rieke, G., et al., 2004, ApJS, 154, 25 

Rodighiero, G., Lari, C, Fadda, D., Franceschini, A., Elbaz, D., 

Cesarsky, C, 2004, A&A, 427, 773 
Scott, S.E., et al., 2002 MNRAS 331, 817 

Scott, S.E., Dunlop, J.S., Serjeant, S., 2006, MNRAS, 370, 1057 
Serjeant, S., et al., 2003a, MNRAS, 344, 887 
Serjeant, S., et al., 2003b, MNRAS, 346, L51 
Serjeant, S., et al., 2004, ApJS 154, 118 

Small, I., Ivison, R.J., Blain, A.W., Kneib, J.-P., 2002, MNRAS 
331, 495 

Small, I, Chapman, S.C., Blain, A.W., Ivison, R.J., 2004, ApJ, 
616, 71 

Tacconi, L.J., et al., 2006, ApJ, 640, 228 
Takagi, T., et al., 2007, MNRAS in press [Paper V] 
van Kampen, E., et al., 2005, MNRAS, 359, 469 
van Kampen, E., et al., 2007, in preparation [Paper X] 
Wagg, J., et al., 2007, MNRAS, 375, 745 
Wang, W.-H., Cowie, L.L., Barger, A.J., 2006, ApJ, 647, 74 
Webb, T.M.A., Brodwin, M., Eales, S., Lilly, S.J., 2004, ApJ 605, 
645 



